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1. Introduction 

At the moment cosmology is one of the fastest changing fields in physics. This fact 
might, on the one hand, be ascribed to the vast amount of new observational data (cf 
fl5| , |16| , |PT1 e.g.), on the other hand there are still fundamental open questions within 
what is nowadays called the cosmological standard model |], p|: Where does the inflaton 
field come from? Is there a something like the cosmological constant A which contributes 
to the dark energy etc.? 

From a theoretical viewpoint one might divide efforts today within cosmology into 
two broad subclasses. Firstly, we have models which extend the standard model to a 
certain amount, inflation pL9|| , e.g., can be viewed as an add-on for the classical 
Friedman-Lemaitre-Robertson- Walker (FLRW) model. All of these models have in 
common that they do not affect the structure of spacetime itself, i.e., they are still 
bound to a four-dimensional Riemannian spacetime and, in addition, do not modify 
the underlying gravity theory, i.e., General Relativity (GR). Secondly, we have models 
which are no longer tied to a four-dimensional Riemannian spacetime and might also 
modify the underlying gravity theory. One example for the latter is the so-called metric- 
affine gauge theory of gravity (MAG), as proposed by Hehl et al in ||. Within this 
gauge theoretical formulation of a gravity theory there are new geometrical quantities, 
torsion T a := Dd a and nonmetricity Q a p := —Dg aj 3, which liberate spacetime from 
its Riemannian structure. In this paper we are going to work within this more general 
framework of a post-Riemannian gravity theory. 

The main reason to consider more general structures within cosmology is the idea 
that new geometrical quantities might shed light on the problems of the cosmological 
standard model, e.g. provide an explanation for the rather artificial introduction of an 
additional scalar field, like the inflaton field. The new quantities couple to the spin, 
shear, and dilation currents of matter, which are supposed to come into play at high 
energy densities, i.e. at early stages of the universe ||, [K|. Another open question, 
which might also be attacked, is the origin of the large amount of dark energy, predicted 
by recent supernova observations. A model which deals with a new kind of dark matter 
interaction within a post-Riemannian theory has been proposed by Tucker and Wang 
in 
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Within this paper, we confine ourselves to a Weyl-Cartan spacetime. This type of 
spacetime can be viewed as a special case of the more general metric-affine framework, 
in which the tracefree part of the nonmetricity Q* a p vanishes by definition. The reason 
to consider this kind of restriction is twofold. Firstly, computations are more feasible 
in a spacetime which is not endowed with the full MAG symmetries. Secondly, the 
Weyl-Cartan spacetime is, unlike the Riemann-Cartan or the Weyl spacetime, still 
able to carry both of the new field strengths nonmetricity and torsion. Let us note 
that the metric-affine framework incorporates all of the above mentioned types of 
spacetimes. The Einstein-Cartan theory, which is formulated in a Riemann-Cartan 
spacetime, represents a viable gravity theory with torsion. By switching off all non- 
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Riemannian quantities in MAG we arrive at GR. 

As becomes clear from the title, this paper stands at the beginning of series of 
articles on cosmological models in alternative gravity theories. In this first article we 
derive the field equations for an enhanced gauge Lagrangian and look for solutions of 
these equations. Thereby we extend earlier joint work with Tresguerres |J and lay the 
foundation for the second article in this series, which will deal with the observational 
consequences of this new cosmological model. 

The plan of this paper is as follows. In section |2] we derive the field equations and 
Noether identities for a new gauge Lagrangian on a formal level. After that we make 
use of computer algebra and provide the explicit form of the field equations and Noether 
identities for a rather general choice of the Weyl 1-form in section In section ^ we 
restrict our considerations to a special case, which leads to a more manageable set of 
field equations. In section [| we look for exact solutions of these equations. We will 
draw our conclusion in section ^ and present some plots for one branch of our model. In 
Appendix A| and [Appendix B| we provide a short introduction into MAG and the Weyl- 



Cartan spacetime. Additionally, we show in [Appendix B| how our model fits into the 



general framework of MAG as proposed in j5| . In [Appendix D| , we provide an overview 
over the units used throughout the preceding sections. Note that we make extensive 
use of differential forms within this paper. A short compilation of symbols used within 



the paper can be found in |Appendix C[ for a more rigorous treatment the reader should 



consult Appendix A of 0]. 

2. Lagrangian, gauge, and matter currents 

In H we considered the following gauge Lagrangian 



V old = ^R a P A V* + Va,%/A *R p a + bZ aP A *RP a . (1) 

i=i 

Where R a p = W a p + Z Q/3 = antisymmetric + symmetric part of the curvature, and 
Va/3 '■= * ($a A ftp). Numbers in parentheses in front of quantities correspond to the 
irreducible decompositions performed in ||. 

Since we are interested in more general gauge Lagrangians (cf equation ( |A.11| ) for 
a very general one proposed in MAG), we are going to extend (|l|). We will perform 
our calculations on the basis of a Weyl-Cartan spacetime, i.e. a spacetime in which the 
tracefree part of the nonmetricity Q a p vanishes, a short introduction into Weyl-Cartan 
spacetime is given in Appendix B . The most obvious extension of (JT|) is given by 



i 



^ = ^ C/ «Q a/3 A^ a . (2) 

7=1 

Since in a Weyl-Cartan space the nonmetricity is reduced to its trace part, i.e. 
Qap = 9apQ = i5 , o / 3 ( 5 7 7 = see equation (|B.1| ), equation (§) now reads 

Vl = c 4 ^Q a/3 A * WQ? a =: cQ aP A (3) 
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Our new Lagrangian now reads 

V = Km + K (4) 
= Einstein-Hilbert + quadratic rotational curvature 

+ quadratic strain curvature + quadratic nonmetricity. (5) 

In contrast to [14]] we included an explicit nonmetricity term in our Lagrangian. Note 
that we have the arbitrary constants Xi a i=i...6, b, c, and the weak gravity coupling 
constant k. The Lagrangian in ([5]) can be viewed as another step towards a better 
understanding of the full MAG Lagrangian as displayed in flA.ll ). Since a treatment 



of the full Lagrangian is computationally not feasible at the moment it is necessary 
to successively study the impact of additional terms in the Lagrangian (a review of 
Lagrangians used in MAG and exact solutions of the corresponding field equations can 
be found in [[Ll|]). Together with the quadratic rotational curvature and quadratic strain 
curvature terms, which were already included in our previous work |J, we now have 
an additional post-Riemannian piece in form of a quadratic nonmetricity term which 
enhances the usual Einstein-Hilbert Lagrangian commonly used in general relativistic 
cosmological models. Note that our Lagrangian in (Q) does not include a term with the 
usual cosmological constant. As we will show in section |] our ansatz in (]|) gives rise to 
an additional constant which, on the level of the field equations, will play the same role as 
the cosmological constant in the standard model. Hence we omit an explicit cosmological 
constant term at this stage. From (pp we can derive the gauge field excitations. They 
read 

M a/3 = - Ac*Q Pa = -c* {g Pa Q\) , (6) 
H a = 0, (7) 

H% = ~ 2^" " 2 S ai HI)W ^ ~ \^ (8) 
The canonical gauge energy-momentum is given by 

E a = e a \V + (e a \ RfP) A + (e a \ T^j A Hp + i (e a \ Q M ) A M Pl 

= e a \V + (eJVO A H /3 1 + ^ (e a \Q^) A (9) 
In contrast to @, we now have a non- vanishing gauge hypermomentum || 
E a j3 = - A Hp - g 0J M^ 

Z4c gfh *Q?>«= cg M *{g^Q\). (10) 
The field equations now turn into 

— E a = S a , (11) 

dH a a - E a a = A, (12) 

g-ylaDH^fl — E[ a/3 ] = T a p. (13) 

| Additional assumptions are marked with an " ^4" . Note that we mark new relations, new with respect 
to 0, with an "AT". 
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Note that in eqs. fll2| ) and (0) we decomposed the second field equation into its trace 
and antisymmetric part (cf |Appendix B| ). Since we are interested in the behaviour 
induced by the new part in (|5|), we will confine ourselves to a non-massive medium 
without spin, i.e. r a/3 = 0. Thus, eq. (|13|) turns into 

g^aDH^] - E[ a p] = 0. (14) 

Because we have not specified a matter Lagrangian, we have to take into account the 
Noether identities (cf |Appendix A| and |Appendix Bp , i.e. 

DY, a = (e a \T?) A E^ - l - (e a \Q) a% + {e a \R m ) A 

+ \{e a \R) A A, (15) 

o" Q /3 = ^rg a pdA + tf( a ASjj), (16) 
=0 [a AEfl. (17) 
We can rewrite eq. (|T6| ) by using (|i~2"D 

= —^9apdE 1 1 + #(a A E/j). (18) 

Note that (p76|)-(p~7|) represent the decomposed second Noether identity in case of a 
vanishing spin current. With fll7|) , eq. (p~5|) turns into 

DY, a ± (e a j T' 3 ) A E^ - i (e a J Q) & A E^ + \ (e a J Q) ^ 

Z o 

+ ^(e a J#)AA. (19) 

Thus, we have to solve (pTTj) , (|12|), (|14]) , and (|T7| ) -([T9| ) in order to obtain a solution for 
our model proposed in (|5|). Now let us investigate the matter sources of our model. For 
a vanishing spin current, the hypermomentum /\ a p becomes proportional to its trace 
part, i.e. the dilation current cf ( |B.4p 

A a/3 = j&tf A\. (20) 

The trace part of the second field equation (|l^) yields 

Kp = \g a p{dH\-E\). (21) 

In contrast to , the dilation current is no longer a conserved quantity since dA a/3 ^ 0. 
From (H|), and (|T0|) we can infer that 

A Q/3 = ~ 9aP {2bd*R\ + cg ltl *{g^Q%)) . (22) 
Because of ( |B.10|) , i.e. i? 7 7 ~ dQ, this equation turns into 

A a/3 = ~ 9aP {bd*dQ\ + cg 7fl * {g^Q v v )) . (23) 

Thus, for our ansatz the hypermomentum A a ^ depends on the nonmetricity and vice 
versa. Note that the second term in (|23| ) depends on the coupling constant introduced 
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in eq. (^j). Now let us specify the remaining quantities in our model. Equation (O) 
forces the components of energy-momentum 3-form to be symmetric, thus we choose|| 

S a = T, afi rf , with S a/3 = diag (fJ,{t),p r (t),Pt(t),Pt(t)) . (24) 

Subsequently, we can calculate the metric stress-energy a a p from eq. (|T8| ) 



o-afs = -\g a p d (c * {g^Q\)) + $ (a A (S /3)7 ?f ) . (25) 

Again we obtained a quantity which depends on the Weyl 1-form, i.e. the trace of 
the nonmetricity. Since we want to compare our model with the cosmological standard 
model, we take the Robertson- Walker line element as starting point of our considerations 

= dt, 1 = - 22 dr, d 2 = S(t) r d6, # = S(t) r sin 6 d<j>, (26) 
v 1 — At 2 

with 

ds 2 = tf 6 <g> # - ^ ® tf 1 - # ® - # ® (27) 

As usual, S'(t) denotes the cosmic scale factor and k = —1,0, or 1 determines whether 
the three-dimensional spatial sections of spacetime are of constant negative, vanishing, 
or positive Riemannian curvature. Following the model proposed in 0, we will choose 
the torsion to be proportional to its vector piece T a ~ ^T a and relate it to the Weyl 
1-form as follows 

T a = -Q A § a . (28) 

The only thing missing for setting up the field equations is a proper ansatz for the 
Weyl 1-form Q. In we were able to derive Q from an ansatz for the potential 
of the hypermomentum A, the so called polarization 2-form P. Here we will adopt 
a slightly different point of view. Since we are interested in the impact of different 
choices of the non-Riemannian quantity Q on cosmology, we will directly prescribe it 
in the following. Besides of the fact that we gain direct control of the post-Riemannian 
features of our model, we circumvent the question which type of matter might generate 
the corresponding hypermomentum. This question and the investigation of models with 
a more sophisticated matter model, like the hyperfTuid of Obukhov et al p |, will be 



postponed to later articles. Let us note that our ansatz in equation (|24D is in general not 
compatible with the energy- momentum obtained in ( ||13|| , eq. (3.28)). Both quantities 
are only equal in special cases like the one we will investigate in section |j. Since we do 
not prescribe a matter Lagrangian and use the Noether identities as constraints on the 
matter variables, our approach could be termed phenomenological as suggested in the 



first part of p3|. 



Here we made use of rj a :— *i) a . 
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3. Field equations and Noether identities 



In this section we will derive the field equations and Noether identities resulting from 
specific choices of the 1-form Q which controls nearly every feature of our model. We 
start with a rather general form of Q, namely 



W S(t) ' 



(29) 



where r) denotes an arbitrary function^ of the radial and the time coordinate, and 
S(t) represents the cosmic scale factor of (0). With the help of computer algebra we 
find that the field equations (|TT|), ([12D , and (|T4]) yield a set of four equations. In order to 
compare these new field equations with the ones derived in || (cf eqs. (40)- (43) therein) 
we write them as follows: 



— (a 4 + a 6 ) k 



+ (a 4 + a 6 ) k 





-1 2\ 



+ 



k 



(30) 




k 

IP 



(|) 2 -^(1-^)|) 



J0+(a 4 + a 6 )J(! 



(|) 2 + & (1-^)|) 



(31) 



Jfe 



(32) 
(33) 



Since we have not specified a matter Lagrangian, we have to be careful with the Noether 
identities (|T3p-(|I7p. Our ansatz (^) yields two equations 

SS 3 (3/i + p r + 2pt) + fiS* - 16ScS£ 2 + 2^ r b (l - kr 2 ) - 8£ f cS 2 £ = 0, (34) 
trrir br (l - kr 2 ) + g. b (2 - 3At 2 ) + 4£, r crS 2 £ + S 4 (ft - p r ) = 0. (35) 

Note that in eq. (|33"D we assumed that a 4 7^ — a^. As one realizes immediately, eqs. 
(^0|)- (|33|) are very similar to the ones we obtained in (||, (40)-(43)). There is only a 
change on the rhs, i.e. the matter side, of the above equations in form of additional 

|| Note that this function is not identical with the one used in eq. (24)). It has a slightly different 
meaning since we use it here directly in our ansatz for the nonmetricity. 



A cosmological model in Weyl-Cartan spacetime I 



8 



terms contributing to the pressure and energy density. As we can see from eqs. (pOj)- 
(|32|), the terms proportional to £ r vanish if we make the same assumptions as in 0, 
i.e. £(t,r) — > £(£). Apart from this feature, there is another, more subtle change in 
(|3T)|)-(|3"3|), i.e. a term of the order £ 2 controlled by the new coupling constant c (cf eq. 
(@))d The Noether identities fl34]) and fl35| ) can be transformed to 

1 c/S" 4 

(/j,-p r - 2pt) , (36) 



2crd^ b_ 



lh- kr 2 ) S 

2 v ; ar 



(2 - 3fcr 2 ) ^ 



(37) 



Comparison of (|37|) with (||, eq. (37)) yields a more sophisticated relation between the 
radial and tangential stresses. 

Let us now extract some more information from the field equations. Addition of 
(0) and (H) yields 



2X| |+(^ 
Subtracting (|32|) from 




2b (1 - kr 2 ) | 



(3? 



5 

+ 2/s (a 4 + a 6 



Let us now combine 
(13 



46 (1 - fcr 2 ) | 



(39) 



(H) « d , , , 



The trace of the energy-momentum reads 
S 7 7 = /j,-p r -2p t 




26(1 -^ 2 ) | 



(40) 



2cr d£ 2 6 /V , 9X d£ 2 . , _\ 



(41) 



Since we encountered a system of coupled PDEs, we will confine us to a special case 
in the following in which the field equations turn into a set of coupled ODEs. At this 
point we would like to note that the above situation is reminiscent to the extensions of 

% Formerly the term of the order £ 2 was controlled by the coupling constant b (cf eq. (Q)), and (0, 
(41)-(43)). 
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the classical FLRW model to anisotropic and inhomogeneous metrical structures. For 
completeness we list the surviving curvature pieces for the ansatz in equation (p9|) 



(42) 
(43) 





SS-S 2 
2S 2 






SS + S 2 
2S 2 


+ Va^, 


^00 — 


- 

"n 


- _(4) V„ _ 
^22 _ 



2S 2 



4. Special case £(t,r) — > ((i) 

In this section we will investigate the interesting special case in which Q, cf eq. 
given by a closed 1-form, i.e. 



Q 



S(t) 



The field equations are now given by 

2 



0,4 + 0,5) k 





5 
5 



A; 

S 2 " 




5/ + ^ 



■^2 I + (a 4 + a 6 ) « 




SI + S 2 



d S ( S\ k . 
Jt[s + (s) + S5 l=a 



is 



(45) 



(46) 



(47) 



(48) 
(49) 



Thus, the function £ contributes to the energy density and pressure in a similar way 
as the function £ in ||. We note that there is no additional contribution from the 
strain curvature in eqs. (|46| ) -(|33| ) , i.e. no term controlled by the coupling constant b 
of our Lagrangian (cf eq. ([]])). This behaviour is explained by the fact that the strain 
curvature vanishes identically for closed 1-forms, like Q from eq. (ffH), in a Weyl-Cartan 
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spacetime. The Noether identities now read: 

p r -pt = 0. 



(50) 
(51) 



In contrast to (0), eq. ( |5T] ) forces the radial stress to be equal to the tangential stress. 
Addition of (|§ and (g§), i.e. eq. ©, yields 



S 



S 



5* 



5 



S 2 



^ [ /x - 3p r + 8c f - 



(52) 



Subtracting (|48| ) from (^) (cf eq. ^) yields 

2 

s* ... 

2%- + 2k (a 4 + a 6 




S 



o I A* + 3 Pr - 16c ( ^ 
Combination of fl5"2"D with the field equations (cf eq. 

2 N 







7737 A 4 — 3p r + 8c — 



3 rft 



c 



5 



7^ 



leads to 



(53) 



k d 
3d* 



~ | S\ + 8C( ^ 



S 7 7 + 8c ( ^ 



const =: H 



(54) 



Thus, we obtained a conserved quantity similar to the one in (0, eq. (47)). The first 
Noether identity (0) takes the form 



d 
dt 



5 4 |/x-8c(| 



^ 1^ 

dt 4 dt 7 



4 (h -//--- 



5 



4 



A = 8c ( - 



2 H 



(55) 



Before we proceed with the search for explicit solutions, we will collect the remaining 
field equations 

2\ 




+ 



k 



\ ( A + - ) -i i<n -(/ (i )K 



(a 4 + a 6 ) k 



c 




5 
5 



+ 



52 



(56) 



S 
S 



-1 2\ 



5 
S 



+ 
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Table 1. Assumptions made up to this point. 



Ansatz/ Assumption Resulting quantity/equation 



Equation 



Tq/3 = 



2* 



— S(t) 17 

a 4 ^ -ae 
E = = 
A 



Affects the form of the second field equation 

Affects the form of the connection 

Controls non-Riemannian features/ 

Affects the form of the field equations 

Controls non-Riemannian features/ 

Simplifies field equations 

Affects the form of the second field equation 

Relation between [i and p r 

Affects the form of the field equations 





(57) 



(5f 



(59) 

^ <s J 

Note that the new constant A in ( |5"7| ) is defined via (p5|). Comparison of ( |5*5| ) with the 
Friedman equation in standard cosmology reveals that A plays the same role as the 
usual cosmological constant. Since we did not include this additional constant in our 
Lagrangian right from the beginning A might be termed induced cosmological constant. 
Now let us exploit the fact that we are allowed to set the constant 5 = 0, which leads 
to an additional constraint, i.e. 

H = 3p r - 8c ( 41 • (60) 



Subsequently eq. (|56f) turns into 



— (a 4 + a 6 ) k 





+ 



S 2 



and the second Noether identity (|59f) now reads 

S ' ^ 2 

4-/i + /^ = 

Note that we collected all assumptions made up to this point in table |]. 




(61) 



(62) 
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5. Solutions 



5.1. A t^O solutions 

We are now going to solve eq. (|58|) for nonvanishing A. Note that eq. (|58|) does not 
depend on the relation between the energy density and pressure and therefore can be 
solved independently. This ODE, after a substitution, turns into a Bernoulli ODE, 
which in turn can be transformed into a linear equation. After this procedure we obtain 
two branches for the scale factor. They read as follows: 

1 



S 



2At 



Ike 



/2At 



2Ax x e 



2V2At 



2Ax 2 



(63) 



where xi, and x 2 are constants. This solution for the scale factor is valid for all three 
possible choices of k. Let us now proceed by fixing the equation of state. 

We will start with the most simple ansatz, i.e. with the introduction of an additional 
constant w into the equation of state, which parameterizes the ratio of the energy density 
and the pressure in our model, 



w fi(t) =p r (t). 

Now let us derive the impact of 
Equation flBTiD yields 



on our set of field equations given by 



ic 



I' 



1 — 3w 

The field equations are now given by 

VI I — I • 4r I 0' "<;) 



(64) 

HUD- 

(65) 




(66) 



3w-l 

Equation (|69D has two non-trivial solution, namely 



c 



with l = const, 



and 



w 



(67) 
(68) 
(69) 

(70) 
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Table 2. Ansatz A ^ 0,w = const. 



Additional constraints 



S from eq. ( p3|) {04 = ag = c = \ = 0}, 
{04 = ag = t = 0, w = 
{a 4 = a 6 = x = 0, ui = 1}, 
{04 = ae = c = 0, w = i}, 
{04 = — ct6, A = 0, w — const} 



Table 3. Ansatz A^0,w = iu(i). 



c 


5 


Additional constraints 


C arbitrary, w = f 2g2j ^_ 3 


5 from eq. (|63 


j) {a 4 = -as, A = 0}, 

{a 4 = -a 6 ,x = 0,w= 1}, 
{a 4 = -ae, 1 = X = 0} 



Solving the remaining field equations with respect to the first solution in eq. (|70|) , 
we obtain constraints among the coupling constants which are summarized in table ^| 
(note that every set of parameters on the rhs corresponds to a solution of the field 
equations). These solutions are not very satisfactory since they either lead to vanishing 
post-Riemannian quantities or, in case of x = c = 0, to a restriction on the Lagrangian 
level. 

Let us switch to another ansatz for the equation of state, namely 

w(t)n(t)= Pr (t). (71) 

Thus, we introduced an additional function into the equation of state which controls the 
relation between the energy density and stresses in a dynamical way. The field equation 
which changes with respect to the set (|66|) -(|69|) , besides of the fact that w is no longer 
a constant, is the Noether identity in eq. (|69|) , which now reads 



(3w - l) 2 



(I) 2 ( I (1 + 3l ° 2 - 4w ) + ^ (! + 3 ^ 2 - H +A=0. (72) 



In case of an arbitrary choice of £, this equation is solved by 

S 2 C 2 K 3 -1 , . 

W = ^2*3-3" (73) 

Reinsertion of this solution for w into the remaining field equations yields additional 
parameter constraints which are summarized in table |3|. As one realizes immediately, 
none of the solutions collected in table |3] is of use for us, since they all lead to unrealistic 
or forbidden restrictions among the coupling constants in our model. Therefore, in the 
following section, we will switch to the case in which the induced cosmological constant 
A vanishes. 
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Table 4. Ansatz A = 0, w 


= tu(t). 




C 


5 


Additional constraints 


C = jj, w = const 


fc ^ 0, 5 from eq. (fr|) 

fc = 0, 5 = const cf cq. @ 

k = 0, S from eq. @ 


A 2 , 

4cz, 2 k+3>i:ix 
10 = 1 

_ 4ct 2 K+^ 2 X 
~~ 4a 2 K+3> f 2 x 


w = c 2 s 2 S-3' ^ arbitrar y 


fc ^ 0, 5 from eq. @ 

fc = 0, 5 = const cf eq. @ 

fc = 0, 5 from eq. @ 


c = 

a ^< 2 x 



5.^. A = solutions 



Solving equation fl5q ) for vanishing A, yields a solution for the scale factor which depends 
on the value of the constant k: 



k ^ : S = ±\/ — {h\ — k 2 (x 2 + t) 2 ) , with x\, x 2 = const, (74) 



= : S = K\Oi S = ± y 2xi (t + x 2 ) , with Xi, x 2 = const. (75) 

Motivated by the results in the previous section for the A 7^ case, we will directly start 
with the more general equation of state as given in ([TI]). The field equations are now 
given by eqs. (|66|)-( |69| ) but with A = 0. The parameter constraints for this solution are 
summarized in the second part of table |j. 

Additionally, we investigated the case in which made use of the old solution for £, 
i.e. as given in ([70]). With this ansatz for £ the second Noether identity, as given in eq. 
(f72f), turns into: 

w = 0. (76) 

Thus, w has to be a constant which, subsequently, can be determined from the remaining 
field equations after choosing the branch for S from eqs. (ffi)-(|f^). The additional 



constraints for this parameter choice are listed in the first part of table [|. Most 
interestingly it turns out that the parameter w, which controls the equation of state, is 
restricted by the choice of a certain set of constants in our theory (cf rhs in table [|). 



6. Conclusion 



In the the last section we have shown that it is possible to find exact solutions of the field 
equations within our model. We were able to generate a rather broad class of solutions 
which allows for a flexible equation of state. We collected the resulting constraints of 
the parameters in our model in tables §-|4|. There seem to be no reasonable solutions in 
case of a non-vanishing induced cosmological constant A, unless one wants to introduce 
strong restrictions on the Lagrangian level. Thus, we are going to focus on the solutions 
with vanishing A in the following. 
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In figure [I] we plotted the scale factor for all three possible values of k and for 
different values of the parameter K\. As becomes clear from the plot at bottom right, 
we have three qualitatively different behaviours depending on the value of k. As in the 
Friedman case the collapsing scenario corresponds to a universe with positive spatial 
curvature. In figure |] we plotted the function Q for the ansatz mentioned in equation 
(fTOf). As stated before Q, the Weyl 1-form controls the non-Riemannian features of our 
model. From the plots it becomes clear that it is possible to construct models in which 
Q vanishes at later times. Thus, the non-Riemmanian quantities die out with time. This 
is a rather desirable behaviour, since the spacetime we are living in nowadays seems to 
be a Riemannian one. At least all experiments carried out so far point into this direction 
0. Nevertheless our model is flexible enough to cope with both situations, i.e. if there 
is evidence for non-Riemannian structures at the present time, we are able to implement 
this fact by modifying our ansatz in ( p9|) and ([45]), respectively. 

In comparison with the usual FLRW model of cosmology we still have three distinct 
cases for the evolution of the scale factor, which correspond to the three different choices 
for k in the ansatz for the metric in equation ( f26| ) . Since one of our field equations (|~)8|) 
is very similar to the Friedman equation in standard cosmology we obtain a similar root 
type behaviour for the scale factor as displayed in figure [I]. As shown in (^) an induced 
cosmological constant leads to inflationary like solutions. In contrast to our old model [H] 
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we were not able to find meaningful parameter constraints for this branch of the model 
(cf . tables |2] and |^) . This drawback might be relaxed in the future if we switch to another 
ansatz for the Weyl 1-form Q. Most interestingly the non-Riemannian quantities lead 
to a contribution to the total energy density of the universe as shown in fl65l). Thus, 
the energy density is no longer a quantity which is determined by the evolution of the 
scale factor only, like in the FLRW scenario. As we will show in the next article of this 
series this contribution might be used to define a new energy density parameter which 
adds to the total energy budget of the universe. Thereby leading to an interesting new 
source for a possible dark energy component. Since the field equations differ from the 
Friedman equations one can expect several observational changes with respect to the 
standard FLRW model. Note that an ansatz with a position dependent Weyl 1-form at 
very early stages of the universe might contribute to the observed inhomogeneities in the 
cosmic microwave background. Although speculative at this time, small inhomogeneities 
in the new geometric quantities might also have served as seeds for structure formation 
at early stages, thereby yielding an interesting supplement to the quantities within the 
standard paradigm. 

Let us summarize that the solutions found above contribute to the collection of 
known exact solutions in MAG, see |TTJ. Additionally, we managed to extend the model 
proposed in M. We provided the foundation for upcoming articles which will deal with 
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the observational consequences of this new model. The most pressing task will be to look 
for realistic parameter choices in order to determine whether the model is in agreement 
with recent observational data. Within the next article of this series [7), we will use the 



supernova data of Perlmutter et al [|T5| and Schmidt et al |16| in order to constrain the 
free parameters in our model. 

Appendix A. MAG in general 

In MAG we have the metric g a/ 3, the coframe d a , and the connection 1-form 
(with values in the Lie algebra of the four- dimensional linear group GL(4,R)) as 
new independent field variables. Here a, (3, . . . = 0, 1,2,3 denote (anholonomic) frame 
indices. Spacetime is described by a metric-affine geometry with the gravitational 
field strengths nonmetricity Q a p := —Dg a/ 3, torsion T a := D$ a , and curvature 
R a P := dT a P — Ta 1 A T^. A Lagrangian formalism for a matter field \l/ minimally 
coupled to the gravitational potentials g a p, $ a , has been set up in ||. The dynamics 
of an ordinary MAG theory is specified by a total Lagrangian 

L = VW^^Q^^iC 3 ) + L mat (^,tf a , (A.i) 

The variation of the action with respect to the independent gauge potentials leads to 
the field equations: 

=0, (A.2) 

DM a/3 - m aP = a aP , (A.3) 
DH a -E a = £«, (A.4) 
DH« P - E°p = AV (A.5) 



Equations ( |A.3| ) and ( |A.4j ) are the generalized Einstein equations with the symmetric 



energy-momentum 4-form a af3 and the canonical energy-momentum 3-form S Q as 
sources. Equation (|A.5 ) is an additional field equation which takes into account 



other aspects of matter, such as spin, shear and dilation currents, represented by the 
hypermomentum A a p. We made use of the definitions of the gauge field excitations, 

rj 9V M AG jja 9V M AG nja f3 9V M AG /. „n 

of the canonical energy-momentum, the metric stress-energy, and the hypermomentum 
current of the gauge fields, 

^:=^f, m<*:=2^, E% = -<P A H p - g M M«\ (A.7) 

and of the canonical energy-momentum, the metric stress-energy, and the 
hypermomentum currents of the matter fields, 

ov a dg a3 6T a P 
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Provided the matter equations QA..2 ) are fulfilled, the following Noether identities hold: 



1 



DT, a = {e a \T?) AE r - (e a \Q 01 ) + (e a J V) A A* 7 , 



(A.9) 

DA a p = g^o - 1 — #* A E^. (A.10) 

They show that the field equation ( |A.3| ) is redundant, thus we only need to take into 
account ( |A.4| ) and ( |A.5| ). 

As suggested in [JTlj], the most general parity conserving quadratic Lagrangian 
expressed in terms of the irreducible pieces of the nonmetricity Q a/ 3, torsion T a , and 
curvature R a p reads 

3 



MAG 



1 r 

2k 



- a R a/3 A nap - 2X V + T a A * ^ aj {I) T a 



J=i 



+ &5 ( {3) Q a7 A^) A *( (4) Q^A^) 



J=2 



-. 6 5 

— R a ^ A * ujj + £ «^ + w 7 tf Q A (e 7 J ® W^) 

P i=i i=i 

9 

+ z 6 $, A (e a J ( 2 >Z^) + zrf* A ( e 7J (/_4) ^ 7 /3)l • (A.H) 



1=7 



The constants entering (|A.11|) are the cosmological constant A, the weak and strong 
coupling constant k and and the 28 dimensionless parameters 

ao, ... ,a 3 ,bi, ... ,b 5 ,oi, ■■■ ,C4,wi, ... ,w 7 ,zi, ... , z 9 . (A.12) 

This Lagrangian and the presently known exact solutions in MAG have been reviewed 
in [11|. We note that this Lagrangian incorporates the one used in section ^| eq. as 
can be seen easily by making the following choice for the constants in (|A.ll ): 



A, a x , . . . , a 3 , 6 3 , h, c 2 , . . . , c 4 , w 7 , z ll . . . , z 3 , z 5 ,...,z 9 = 0. (A.13) 

In order to obtain exactly the form of (f|), one has to perform the additional 
substitutions: 



a 



-X, w u . . . , w 6 -> -2pa u -2pa 6 , £> 4 -> c, z 4 -> -2pb. (A. 14) 



In table [Al| we collected some of symbols defined within this appendix. 

+ [A] =lcngth- 2 , [k] =lcngth 2 , [p] = [ft] = [c] = 1. 
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Table Al. Summary of definitions made in Appendix A 



Potentials Field strengths Excitations 



Gauge currents 



9af3 

r 13 



Qa{3 

R 



= -Dg a0 M a @ 
= D§ a H a 
="D"rJ H a fi 



o OV 

OV 
OT" 
OV 



m af3 



2 dV 

dg a p 
OV 
00" 

ov 

0T a 



Appendix B. Weyl-Cartan spacetime 

The Weyl-Cartan spacetime (Y n ) is a special case of the general metric-affine geometry 
in which the tracefree part ^f Q/ g of the nonmetricity Q a p vanishes. Thus, the whole 
nonmetricity is proportional to its trace part, i.e. the Weyl 1-form Q := \Q a a , 

Qa/3 = 9af3 Q = -9a/3 <3V (B.l) 

Therefore the general MAG connection reduces to 

r a/3 = ]zdg aP + (e [a \dgp h ) $ 7 + e [a \Cp] - - (e a Je (3 JC 7 ) tf 7 

- e [a J + - {e a \ ep\ T 7 ) $ 7 + -g a p Q + (e [a \ Q) 



rJJ - ^\ Tfl + 9 (e Q J e p \ r 7 ) V + ^g a p Q + (e [a \ Q) 6 



(B.2) 
(B.3) 

Thus, it does not include any more a symmetric tracefree part. Let us now recall the 
definition of the material hypermomentum A a p given in ( |A.8|) . Due to the absence of a 
symmetric tracefree piece in (|B.3|) , A Q p decomposes as follows 



A a p = antisymmetric piece + trace piece 
= r a(3 + i# Q/3 A = T a(3 + -g a p A 7 7 
= spin current + dilation current. 
According to QB.4|) the second MAG field equation (|A.5|) decomposes into 

dH a a — E a a = A, 
g^aDH 1 $ — E[ a p] = T a p, 

while the first field equation is still given by (|A.4j) . Additionally, we can decompose the 
second Noether identity ( |A.10|) into 



(B.4) 

(B.5) 
(B.6) 



-g a p dA + $ (a A Ep) 



(J, 



a/3, 



(B.7) 

Dr a p + Q A r af3 + [a A Ep] = 0. (B.8) 

Thus, the first Noether identity ( |A.9| ) with inserted Weyl 1-form and hypermomentum 
reads 

DE a =(e a \T^)AEp-^(e a \Q)a%+(e a \R [ p l] )Ar^ + ^(e a \R)AA.(B.9) 
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Table CI. Operators. 



Operation 



Symbol Input 



p-form A g-form 
vector J p-form 
d p-form 
*p-form 



Output 



Exterior multiplication 
Interior multiplication 
Exterior derivative 
Hodge star in a n-dimen. space 



(p + q)-form 
(p — l)-form 
(p + l)-form 
(ri — p)-form 



Table Dl. Natural units. 



[energy] 


[mass] 


[time] 


[length] 


length -1 


length -1 


length 


length 



Finally, we note that in a Y n spacetime the symmetric part of the curvature R( a j3) = Z a p , 
i.e. the strain curvature, reduces to the trace part 

Z a p = ^g a pR = -gapR 1 ^ = ^g a p dQ. (B.10) 
Appendix C. Differential geometric formalism 

We assume a connected n-dimensional differential manifold Y n as underlying structure 
throughout the paper. A vector basis of its tangent space T p Y n is denoted by e a , which 
is dual (i.e. e a Ji?^ = (5>f) to the basis ft 01 of the cotangent space T*Y n . A p-form H can 
be expanded with respect to this basis as follows 

S=iH ft ... /%1 ^A...A^. (C.l) 

Table |Cl| provides a rough overview of the operators used throughout the paper. For a 
more comprehensive treatment the reader should consult Q or section 3, and Appendix 
Aof§. 



Appendix D. Units 



In this work we made use of natural units, i.e. h — c—1 (cf table Dl). Additionally, we 



have to be careful with the coupling constants and the coordinates within the coframe. 
In order to keep things as clear as possible, we provide a list of the quantities emerging 
throughout all sections in table D2. Note that [d] = 1 and [*] = length" -2 -", where n = 
dimension of the spacetime, p = degree of the differential form on which * acts. 
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Quantities 


1 


Length 


Gauge potentials 




[B»] 


Gauge field strengths 


[Qaff], [Rap] 


[T a ] 


Gauge field excitations 


[M a %[H a fi ] 


[Hat 1 


Gauge field currents 




[EaV 1 


Matter currents 




Pa]" 1 



Coordinates 
Functions 

Miscellany 
Constants 



W,r)],[C(t)] 



[t] 

[S(t)],Ht)]-i,[p r (t)]-i 

[Pt(t)]j 

[E a/3 ]"4 

[k]3, [A]-i, [ci]~i, [c]-2 



eq 
eq 
eq 
eq 
eq 
eq 



0) 

(75) 
(6c) 
(64) 
(7C) 
(72) 



[*l], [X2[ 
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